Cosmic microwave background bispectrum of tensor passive modes induced from 

primordial magnetic fields 
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If the seed magnetic fields exist in the early Universe, tensor components of their anisotropic 
stresses are not compensated prior to neutrino decoupling and the tensor metric perturbations 
generated from them survive passively. Consequently, due to the decay of these metric perturbations 
after recombination, the so-called integrated Sachs-Wolfe effect, the large-scale fluctuations of CMB 
radiation are significantly boosted. This kind of CMB anisotropy is called the "tensor passive 
mode." Because these fluctuations deviate largely from the Gaussian statistics due to the quadratic 
dependence on the strength of the Gaussian magnetic field, not only the power spectrum but also 
the higher-order correlations have reasonable signals. With these motives, we compute the CMB 
bispectrum induced by this mode. When the magnetic spectrum obeys a nearly scale-invariant 
shape, we obtain an estimation of a typical value of the normalized reduced bispectrum as £i(£i + 
1)^3(^3+1) [6< x ^3 ] ~ (130 — 6) x 10~ 16 (£>iM P c/4.7nG) 6 depending on the energy scale of the magnetic 
field production from 10 14 GeV to 10 3 GeV. Here, BiMpc is the strength of the primordial magnetic 
field smoothed on IMpc. From the above estimation and the current observational constraint 
on the primordial non-Gaussianity, we get a rough constraint on the magnetic field strength as 
BiMpc < 2.6 - 4.4nG. 

PACS numbers: 98.70.Vc, 98.62.En, 98.80.Es 



I. INTRODUCTION 

Cosmological observations have suggested the exis- 
tence of micro-Gauss strength magnetic fields in galaxies 
and clusters of galaxies at the present Universe. As their 
origin, many researchers have discussed the possibility 
of generating the seed fields in the early Universe (e.g. 
[EH). These scenarios have been verified by constrain- 
ing the strength of the primordial magnetic fields (PMFs) 
through the effect on CMB fluctuations. 

Conventional studies have provided upper bounds on 
PMFs with the two point correlations (power spectra) 
of the CMB temperature and polarization anisotropics 
[1, [ij]. On the other hand, taking into account the 
CMB three-point correlations (bispectra), which have a 
nonzero value because the CMB fluctuations are sourced 
from the quadratic (non-Gaussian) terms of the stochas- 
tic (Gaussian) PMFs, some new consequences have been 
obtained. In Refs. [5-7], the authors evaluated the con- 
tribution of the scalar modes at large scale with sev- 
eral approximations, such as the Sachs- Wolfe limit, and 
roughly estimated the upper limit on the PMF strength. 
In our previous papers [H, we computed the effect 
of the vector modes without neglecting the complicated 
angular dependence, and obtained tighter bounds due 
to the dominant contribution at small scale induced by 
the Doppler and the integrated Sachs- Wolfe (ISW) effects 
[Tol lllj . However, if the gravitational waves are gener- 
ated from the PMF anisotropic stresses uncompensated 
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prior to neutrino decoupling, these superhorizon modes 
survive passively and the decay of their modes after re- 
combination amplifies the CMB anisotropics through the 
ISW effect [l^ . This type of fluctuation is called the 
"tensor passive mode" and it is expected that the CMB 
bispectrum of this mode has the most dominant signal 
at large scales, as inferred from the power spectrum [13j . 
Therefore, in this paper, we investigate the exact CMB 
bispectrum of tensor passive modes induced from PMFs 
and place a new constraint on the strength of PMFs. In 
the calculation, because there are complicated angular 
integrals as there are in the vector mode case, we apply 
our computation approach, as discussed in Ref. [|j]. 

This paper is organized as follows. In the next sec- 
tion, we formulate the CMB bispectrum of tensor pas- 
sive modes induced from PMFs. In Sec. IIII1 we show 
our result for the CMB bispectrum and the limit on the 
strength of PMFs, and give a discussion. 



II. FORMULATION OF TENSOR 
BISPECTRUM INDUCED FROM PMFS 

Let us consider the stochastic PMFs B b (x, r) on 
the Friedmann- Robertson- Walker and small perturbativc 
metric as 

ds 2 = a(T) 2 [-dT 2 + 2h ob dTdx b + {S bc + h bc )dx b dx c ] . (1) 

Here a is a scale factor and r is a conformal time. In this 
space-time, the PMF evolves as B b (x, r) = B b (x)/a 2 . 
Then the spatial components of the PMF's energy mo- 
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mentum tensors are given by 
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B 2 (x) 



S b c -B b (it)B c (x) 
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p 7 (A B <5 b c + U b Bc 



where B 2 = B b B b and we use the photon energy density 
p 7 (oc a -4 ) for normalization. In the following discussion, 
the index is lowered by 6 bc , and the summation is implied 
for repeated indices. 



A. Bispectrum of the tensor anisotropic stress 
fluctuations 



less polarization vector satisfying the orthogonal condi- 
tion, k a e { a ±1] = 0, and cr(= ±1) expresses the helicity of 
the polarization vector. In general, the magnetic power 
spectrum should contain an asymmetric helical term [l4l - 
llq |. However, we assume the magnetic fields are isotropic 
and homogeneous, for simplicity; hence, this effect is ne- 
glected in Eq. (|4]). Because the production mechanism 
of PMFs remains to be done, we use a simple power-law 
form as the power spectrum: 
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The Fourier component of Usab is given by the convo- 
lution of the PMFs as 

1 f r/ 3 k' 

lWk) = -- Ba (k')B b (k-k') , (3) 

47rp 7 ,o J (2tt) j 

where p 7i a = p n a denotes the present energy density 
of photons. If B a (x) obeys the Gaussian statistics, the 
power spectrum of the PMFs Pe(fc) is defined by 

(B a (k)B b (p)) = (2nf^^P ab (k)5{k + p) , (4) 
with a projection tensor 

P ab (k) = = ^ ~ hh , (5) 

<r=±l 



where -BiMpc denotes the magnetic field strength 
smoothed on a scale IMpc, fciMpc = 27rMpc _1 , and ub 
is a spectral index. 

With a transverse and traceless polarization tensor 
[l7| , e^ 2 ^ = y/2ea the anisotropic stress fluctua- 

tion is decomposed into two helicity states of the tensor 
mode as 



lWk) = E n st(k)C(k) 



A=±2 

which is inversely converted into 



ng 2) (k) = - e ir ; (k)n Bab (k) . 



(T2), 
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which comes from the divergenceless of the PMF. Here k 
denotes a unit vector, is a normalized divergence- 



(±2) 

From the above equations, the bispectrum of is 
symmetrically formed as 
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where A means two helicities: Ai,A2,A3 = ±2, and the 
curly brackets denote the symmetric 7 terms under the 
permutations of indices: a f> 6, c f> d, or e f> /. Be- 
cause the anisotropic stress fluctuation depends quadrati- 
cally on the Gaussian magnetic fields as shown in Eq. (|3|) , 
the statistics of their tensor modes given by © is highly 
non-Gaussian. Hence, the bispectrum of Eq. ([§]) also has 
a nonzero value and induces the finite CMB bispectrum. 



B. CMB temperature bispectrum of tensor passive 
modes 



As is well known, the gravitational potential of ten- 
sor modes can be generated from anisotropic stresses via 
the Einstein equation. If PMFs exist, the anisotropic 
stresses, as mentioned in the previous subsection, also 
behave as a source. In general, after neutrino decou- 
pling, the anisotropic stresses of PMFs vanish via the 
compensation of those of neutrinos. However, prior to 
this epoch, there is no compensation process due to the 
absence of the neutrino anisotropic stresses. Hence, from 
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the Einstein equation, we find the evolution equation of the tensor-mode metric perturbations as 
I 



M ±2 '"(k,r)+2-^ ±2 )'(k,r) + fc 2 ^ ±2 )(k,r)« ( ^G^ngf } (k) (r E ,< r < r v ) (1Q) 
a [0 (t >t v ) 

I 



where r v and tb are the conformal times at neutrino 
decoupling and the generation of the PMF, respectively, 
and ' denotes a derivative of conformal time. Here ft/ ±2 ) 
is given by 1 

/ l ( ±2 )(k,T) = ieif ) (k)/i a6 (k,r) . (11) 

From Eq. (fTU)) . we find a super horizon solution of the 
tensor metric perturbation as [H, [H| 

h^(k) » hW(k,T v ) « 6i? 7 ln ng 2) (k) , (12) 

where _R 7 ~ 0.6 is the ratio by the energy density of 
photons to all relativistic particles for t < t„. 

The CMB temperature fluctuation is expanded into 
spherical harmonics as A ^ n ^ = ^ hn ai m Y^ m {vL). The 
ag m sourced from the initial tensor perturbations (|12|) 
can be expressed as [l7| 

at m = {-if f ^%(k) £ h l%W > ( 13 ) 

J 77 A=±2 

h^\k) = J d 2 k/,( ±2 )(k) T2 r; m (k) . (u) 

where 7^(fc) denotes the transfer function of tensor 
modes. Because the solution of the magnetic passive 

1 h( ±2 > is equal to 2H T of Refs. [HQ]!. 



mode (|12p. if any, would dominate the tensor- mode per- 
turbation, the evolution of tensor modes after their cre- 
ation is almost identical to the standard cosmological one 
without anisotropic stress sources. Therefore, we can use 
the standard cosmological tensor-mode transfer function 

QJlol. 

The CMB angle-averaged bispectrum is given by 

rci\ra 2 m 3 ^ ' \n— 1 / 

where the bracket denotes the Wigner-3j symbol. 

In order to calculate the bispectrum of ai m given by 
Eq. (fT3"|) . we rewrite all angular dependencies in Eq. © 
in terms of the spin-weighted spherical harmonics with 
the notation as [l7| 

4 ±1} (*) = 4 T1) *(?) = Tl«:±i^(f) , 

, 4. m (16) 

"a "a = Y " m ' ' 

We then express the angular integrals of the spin spher- 
ical harmonics with the Wigner-3j symbols, and sum up 
these Wigner-3j symbols over the azimuthal quantum 
numbers in the same manner as in Ref. @. Then, we 
obtain the final form of the bispectrum as 
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where is the spherical Bessel function, fc^i is the matrices in the curly brackets denote the Wigner-6j and 

Alfven-wave damping length scale, the 2x3 and 3x3 
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9j symbols, respectively, and 
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As shown in Eq. (|17[) , the bispectrum depends on tb- 
Although the production mechanism of PMFs is unclear 
and still being discussed, we assume that PMFs arise 
sometime between the energy scale of any grand uni- 
fication theory and the electroweak transition. Hence, 
in the computation of the CMB bispectrum, we con- 
sider two corresponding values: t u /tb ~ 10 17 , 10 6 . This 
leads to a factor of 23 difference in the amplitude of the 
CMB bispectrum due to the logarithmic dependence on 
t u /tb- Therefore, due to the sextuplicate dependence of 
the CMB bispectrum on the magnetic strength, there is 
a model-dependent factor 23 1/6 ~ 1.7 in bounds with the 
PMF strength. 



III. NUMERICAL RESULTS AND DISCUSSION 

Following the final expression (fT?)) . we compute the 
CMB temperature bispectrum of tensor passive modes 
numerically 2 . 

In Fig. [TJ we describe the reduced bispectra of temper- 
ature fluctuations induced by the PMFs defined as [2l[ 
b iltata = (illlr'Be^, for 4 = 4 = 4- From 
the red solid lines, we can find that the enhancement 
at £ < 100 due to the ISW effect gives the dominant sig- 
nal like in the angular power spectrum Ce [HI, [22[ . The 

amplitude £ ~ 0(1) is comparable to C^ 2 because the 
power-law suppression of the Wigner symbols like the 
vector mode (8| is not effective at small £. 

In Fig. [2j we also show b^ 2 ^ 3 with respect to 4 for 
£\ = £2- From this figure, for tib = —2.9, the normalized 
reduced bispectrum is evaluated as 



4(4 + 1)4(4 + 1) I h^i 3 
~ (130-6) x 10~ 16 



-BlMpc\ 

4.7nG J 



(18) 



10 



17 



where the factor 130 corresponds to the t u /tb 
case and 6 corresponds to 10 6 . It is also clear that bg^^ 
for us ~ — 3 dominates in 4 = 4 ^ 4> Comparing this 
with the approximate expression of the bispectrum of 



2 Unlike the case of the vector mode bispectrum calculation Q, 
we do not use the thin LSS approximation because the temper- 
ature anisotropies from tensor modes are nonlocal. To check 
our numerical calculation, we computed the CMB power spec- 
trum of the tensor magnetic passive mode using the same method 
described in the main text, namely, by expanding the nonlinear 
convolution of magnetic anisotropic stress with the spin- weighted 
spherical harmonics. We observe that our results are consistent 
with the previous results [T3l . 
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FIG. 1: (color online). Absolute values of the normalized re- 
duced CMB bispectra given by Eq. |17J for h = £2 = h- The 
lines correspond to the spectra of tensor passive modes (red 
solid lines), vector modes || (green dashed line), and primor- 
dial non-Gaussianity with /nl^ = 5 [2l|] (blue dotted line). 
The PMF parameters are fixed to BiMpc = 4.7nG, ns = —2.9, 
and t v /tb = 10 17 (upper line) and 10 6 (lower line), and the 
other cosmological parameters are equal to the mean values 
limited from the Wilkinson Microwave Anisotropy Probe 7-yr 
data SI- 



local-type primordial non-Gaussianity in curvature per- 
turbations as [24[ 

4(4 + 1)4(4 + i)bt lta ts ~ 4 x io- 18 /NL cal , (19) 

the relation between the magnitudes of the PMF and 
the nonlinearity parameter of the local-type configuration 
/^L al is derived as 



-BlMpc 

InG 



(1.22 - 2.04)|/^ al | 



local 1 1/6 



(20) 



Using the above equation, we can obtain the upper 
bound on the PMF strength. As shown in Fig. [IJ because 
the tensor bispectrum is highly damped for £ > 100, we 
should use an upper bound on /^l 3,1 obtained by the 
current observational data for £ < 100, namely /nl 3,1 < 
100 [25[ . This value is consistent with a simple prediction 
from the cosmic variance 21( . From this value, we derive 
-BiMpc < 2.6 — 4.4nG These are 4—2 times stronger than 
vector- modes bounds [|[ . 

In this paper, we study the CMB temperature bispec- 
trum generated from the tensor anisotropic stresses of 
PMFs and find a new constraint on the magnetic field 
magnitude when the PMF spectrum is close to a scale- 
invariant shape. Although there is a touch of uncertainty 
in the production epoch of PMFs, this bound is tighter 
than ones obtained by the CMB power spectra 0, Hj]. 
Although this limit is weaker than a rough bound from 
only the scalar passive modes @ due to the rapid decay 
of the tensor bispectrum at small scales, the significant 
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FIG. 2: (color online). Absolute values of the normalized 
reduced CMB bispectra given by Eq. (|17|l and generated from 
primordial non-Gaussianity in curvature perturbations given 
by Eq. (|19[) as a function of I3 with l\ = £2 . Each parameter 
is identical to the values defined in Fig. [1] 



amplitude at large scales will have a drastic impact on 
the precise calculation of the limit on PMFs, including 
the scalar, vector, and tensor-mode contributions. 

In our previous studies and the above analysis, we 
find that tensor (vector) modes dominate at large (small) 



scale, not only in the power spectrum but also in the bis- 
pectrum. It is also expected that the scalar mode dom- 
inates at the intermediate scale. Therefore, using this 
scale-dependent property, we will also constrain a spec- 
tral index of the PMF spectrum in addition to the mag- 
netic strength. These reasonable bounds will be obtained 
by considering the CMB temperature and polarization 
bispectrum of autocorrelations and cross-correlations be- 
tween scalar, vector, and tensor modes in the estimation 
of the signal-to-noise ratio. 
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